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ABSTRACT 

We study the relationship between the field star formation and cluster formation properties in a 
large sample of nearby dwarf galaxies. We use optical data from the Hubble Space Telescope and from 
ground-based telescopes to derive the ages and masses of the young (i age < 100 Myr) cluster sample. 
Our data provides the first constraints on two proposed relationships between the star formation rate 
of galaxies and the properties of their cluster systems in the low star formation rate regime. The 
data show broad agreement with these relationships, but significant galaxy-to-galaxy scatter exists. 
In part, this scatter can be accounted for by simulating the small number of clusters detected from 
stochastically sampling the cluster mass function. However, this stochasticity does not fully account 
for the observed scatter in our data suggesting there may be true variations in the fraction of stars 
formed in clusters in dwarf galaxies. Comparison of the cluster formation and the brightest cluster in 
our sample galaxies also provide constraints on cluster destruction models. 



1. INTRODUCTION 

In galaxies, the relationship between stars in bound 
clusters and those that populate the field is poorly under- 
stood. The formation of both populations is connected, 
whether stars form in clusters and dissolve to create 
field stars, or field stars form concurrently with clusters. 
Until recently, it was widely accepted th at the major- 
ity of stars form in a clu ster environment (|Lada fe Ladal 
2003; Porras ct al. 2003). However, infrared surveys have 
shown evidence that star formation occurs in a contin- 
uum of stellar d ensities, including formation of unboun d 
field stars (e.g.. lEvans et al.ll2009t iBressert et al.ll2010D . 
We study the behavior of field and clustered star forma- 
tion in dwarf galaxies by quantifying the fraction of stars 
that populate star clusters. 

A quantity that directly probes the fraction of stars lo- 
cated in clusters is the Cluster Formation Efficiency (F), 
defined as the Cluster Formation Rate divided by the 
Star Formation Rate (CFR/SFR; iBastianl [20081) . This 
quantity has been observed to positively correlate with 
the star formation rate density "Esfr (SFR/kpc 2 ), sug- 
gesting that higher star formation rates provide an envi- 
ronme nt beneficial for cluster formation. Goddar d et al.l 
(2010) first found this trend to be a power-law relation- 
ship for spirals and sta rburst dw arf galaxies (s e e also 
lLarsen fc Richtlerl I200CJI) . Later, lAdamo et at] (|2011[ ) 
found that the correlation extended to galaxies with 
extreme star formation (e.g., starb urst nucl e i and blue 
compact galaxies). Howev er, both IBastianl ([2008D and 
ISilva- Villa fc Larsenl ([201 ID found no clear trend in T 
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with the star formation rate in their studies of spiral and 
starbursting galaxies. In this paper we will test whether 
this trend is seen at lower SFR densities and whether 
significant scatter exists around the trend. 

Determination of the star cluster formation efficiency 
in low SFR density galaxies has complications due to 
cluster populations which are dominated by low mass 
clusters. Recent studies have found that estimated clus- 
ter parameters when derived by comparing observed pho- 
tometry to synthetic spectra show a large spread in de- 
rived ages an d masses for clusters below a few times 



10 4 Mq (e.g.. lAnders et al ; ll2004t IPiskunov et al.l 
Fouesneau fc Lancon l I2010t IPopescu fc Hanson! 



2009; 



2010; 



Silva- Villa fc Larsenl 1201 11) 7 The spread in derived ages 



and masses results from the random presence or absence 
of a few post main-sequence stars that dominate the 
light of the cluster. For instance, a young cluster of 
10 4 M Q is expected to have on average one supergiant, 
but could end up with multiple or no supergiants. As 
a result, adding a luminous red supergiant to the aver- 
age prediction will shift the fitted population towards 
older ages. Bayesian methods are now available that 
can explicitly take into account these stochastic fluc- 
tuati ons of integrated color s and fluxes of star clusters 
(Fouesnea u fc Lanconl [2010'). These effects need to be 
taken into account when determining the cluster forma- 
tion rate of low mass systems. 

Another complication of measuring T for our sample 
is due to the inconsistency of different star formation 
rate indicators in dwarf galaxies. A galaxy's SFR can be 
estimated with many different tracers (Ha, integrated 
far-infrared, UV, resolved star Color-Magnitude Dia- 
gram (CMD) fitting, etc.). Unfortunately, the integrated 
flux tracers, although useful, yield SFR approximations 
based on stars of many ages and become in consistent 
with each other in low mass g alaxies (Mcurcr ct al. 2009; 
Lee et al.1 l2009bl: jBoselli eTall [20091 : .Hunter et al.H201C 



Weisz et al.ll2012h . Ideally, knowledge of every individ- 



ual star's age and mass would provide an accurate SFR. 
However, the closest method to this is reconstructing star 



formation histories using individually resolved stars in a 
galaxy. Utilizing HST's high image resolution capabil- 
ities, the star formation histories can be calculated in 
nearby galaxies (D < lOMpc). Impressive temporal res- 
olution at young ages (t agc < 100 Myr) can be obtained 
and can provide an excellent foundat ion against which 
young cluster ages can be compared dWeisz et al.ll2008t 
I'McQuinn etHI l20Tlt iSilva- Villa fc Larsenl 1201111 The 
mass and age distributions of star clusters combined with 
CMD-derived field star formation histories can produce 
a direct measurement of the number of stars that remain 
in bound star clusters, T. 

Despite the the difficulties associated with analyzing 
dwarf galaxy properties, there are many advantages. 
Dwarf galaxies are abundant in the local universe and 
lack spiral structure. Also, low extinction plus low sur- 
face brightness make finding star clusters relatively easy. 
While few clusters are formed in individual low star 
formation rate dwarfs, the large number of local dwarf 
galaxies enables construction of a statistically significant 
cluster sample from many galaxies. These properties 
make dwarfs an ideal laboratory for identifying star clus- 
ters and studying individual resolved stars. 

In this work we will test existing relationships found 
between the cluster and field star populations at young 
(^agc < 100 Myr) ages for a sample of low star formation 
rate density, dwarf galaxies, and we examine in detail if 
the cluster sample is consistent with a universal cluster 
formation mechanism. 

2. DATA 

In this section we present the data and analysis used to 
derive cluster and star formation rates in a large sample 
of dwarf galaxies. We identify clusters using HST data, 
and derive their ages and masses using a combination 
of HST and ground-based photometry. We also use the 
HST imaging of individual stars to derive recent star 
formation histories of each galaxy. 

2.1. Sample & Observations 

The sample consists of 37 dwarf galaxies selected from 
the ACS Nearby Galaxy Surve y Treasury (ANGST) 
project of iDalcanton et all (|2009t) . ANGST is a volume- 
limited survey and and therefore it contains a large num- 
ber of dwarf galaxies. We restricted the ANGST sample 
to galaxies fainter than Mb = —18 in order eliminate 
non-dwarf galaxies. We further restricted the sample to 
contain galaxies that have available Advanced Camera 
for Surveys (ACS) data (some ANGST data was taken 
with WFPC2), and that are observable from the northern 
hemisphere. All of these restrictions produce a sample of 
37 dwarf galaxies from the original 69 of the ANGST par- 
ent sample; properties of the final sample can be found 
in Table 1. The distance range of the sample is ~2-4.5 
Mpc; at these distances the majority of star clusters are 
partially resolved into individual stars at the resolution of 
HST. Crowding issues near the core of the clusters result 
in an incomplete census of the resolved stars and there- 
fore make it impossible to derive accurate cluster ages 
and masses from CMD fitting methods. Therefore, we 
supplemented ANGST data with ground-based data to 
obtain a Johnson UBV and Cousins RI data set for each 
galaxy in our sample that contained clusters, thereby en- 



abling accurate integrated-light age and mass determina- 
tions. 

The majority of galaxies in our sample have ANGST 
ACS Wide Field Camera (WFC) data taken in the 
F555W/F606W (~V) and F814W (~I). We therefore 
obtained ground based data in the U, B and R filters. 
The UBR ground-based data came from both the BOK 
telescope atop Kitt Peak and WIRO. The 2.1 meter BOK 
telescope, operated by Steward Observa tory, in the 90 
prime configuration ([Williams et al.1 12004) has a blue op- 
timized Anti-Reflective (AR) coated CCD which facili- 
tates deep U images with significantly shorter exposure 
times than with traditional non-coated CCDs. The BOK 
data were taken over three runs: November 2009, May 
2010, and December 2010. The remaining B and all of 
the R data were taken at the WIRO 2.3 meter telescope 
in the prime focus configuration over the months of April 
2010, March 2010, September 2010, March 2011, and 
July 2011. Two galaxies, DD078 and UGCA292, have 
ACS WFC data in the F475W and F814W filters. For 
these galaxies, we obtained broadband V images of these 
galaxies at the Wyoming InfraRed Observatory (WIRO) 
in order to complete the UBVRI data set. There are 
also two galaxies whose full data set was not obtained 
(GR8 and UGC5442); these galaxies are not included in 
the analysis. 

2.2. Cluster Identification 

The cluster catalog for the full ANGST sample will be 
presented in a separate work (Johnson et al. 2012; sub- 
mitted); we briefly summarize how clusters were identi- 
fied here. Cluster identification was performed by four 
members of our team via visual inspection of HST color 
images. Potential objects were chosen based on grouping 
of stars, spherical symmetry, colors th at are consistent 
with single stellar p opulation models (M arigo et al .1120081 ; 
IGirardi et alj|20 10), and a lack of structure (spiral arms 
or otherwise) that might indicate a background galaxy. 
However, since a cut based on colors that are inconsistent 
with non-stochastic stellar models could potentially re- 
move young, low-mass clusters from the sample, we chose 
only to cut clusters who's color severely deviate (greater 
than one magnitude) from the model colors. After initial 
candidates were chosen from the images, all candidates 
were re-examined and ranked according to the likelihood 
that they were bound star clusters. Background elliptical 
galaxies are the most likely contaminant in our sample 
as these can be difficult to distinguish from older clus- 
ters if they fall on the face of the dwarf galaxy. However, 
since we only aim to constrain the formation properties 
of young clusters (t age < 100 Myr), the red colors of po- 
tential background galaxies will typically exclude these 
objects from the final sample of young clusters. 

In total, we detect 144 clusters, 44 of which have an 
age less than 100 Myr in our 37 galaxy sample. Con- 
versely, 20 galaxies in our sample contain no clusters at 
any age. We have tabulated the photometric properties 
of the clusters with ages less than 100 Myr in Table 4 lo- 
cated in the appendix. We include this table to facilitate 
comparison of our cluster fitting methods (discussed in 
Section 2.4) to others' methods. 

2.3. Photometry 







TABLE 1 
Dwarf Galaxy Parameters 








Galaxy 


Alternative 


RA 


Dec 


a 


b 


M B 


Dist 


12+log(0/H) 




Name 


(J2000.0) 


(J2000.0) 


(") 


(") 


(mag) 


(Mpc) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


ESO540-G030 


KDG2 


00 49 21.1 


-18 04 28 


84 


74 


-11.29 


3.33 




NGC404 


UGC0718 


01 09 26.9 


35 43 03 


105 


105 


-16.25 


3.05 




KKH37 




06 47 45.8 


80 07 26 


52 


10 


-11.26 


3.26 




NGC2366 


UGC3851.DD042 


07 28 52.0 


69 12 19 


308 


158 


-15.85 


3.21 


7.95 ± 0.05 b 


UGCA133 


DD044 


07 34 11.3 


66 53 10 


108 


76 


-11.89 


3.10 




UGC4305 


HoII,DDO50 


08 19 05.9 


70 42 51 


278 


232 


-16.57 


3.38 


7.72 ± 0.14 a 7.92 ± 0.10 b 


M81dwA 


KDG52 


08 23 56.0 


71 01 46 


39 


39 


-11.37 


3.44 




UGC4459 


DD053 


08 34 06.5 


66 10 45 


67 


56 


-13.23 


3.61 


7.60 ± 0.11 a 7.82 ± 0.09 b 


UGC5139 


HoI,DD063 


09 40 28.2 


71 11 11 


132 


110 


-14.26 


3.90 


7.60 ± 0.11 a 8.00 ± 0.10 b 


[FM2000] 1 




09 45 25.6 


68 45 27 


41 


11 


-10.16 


3.53 




BK3N 




09 53 48.5 


68 58 09 


20 


20 


-9.23 


3.86 




KDG61 




09 57 02.7 


68 35 30 


107 


60 


-12.54 


3.49 


8.35 ± 0.05 b 


Arp's loop 




09 57 29.0 


69 16 20 


68 


68 


-11.16 


3.78 




UGC5336 


HoIX,DD066,KDG62 


09 57 32.4 


69 02 35 


124 


90 


-13.31 


3.61 


8.14 ± 0.11 a 8.65 ± 0.25 b 


LEDA166101 


KK98.77 


09 50 10.0 


67 30 24 


110 


76 


-11.42 


3.55 




UGC5428 


DD071,KDG63 


10 05 07.3 


66 33 18 


98 


81 


-11.71 


3.53 


7.4 C 


UGC5442 


KDG64 


10 07 01.9 


67 49 39 


98 


62 


-12.32 


3.72 


7.4 C 


IKN 




10 08 05.9 


68 23 57 


90 


78 


-10.84 


3.61 


7.3 C 


[HS98] 117 




10 21 25.2 


71 06 58 


106 


64 


-11.51 


3.82 




DD078 




10 26 27.9 


67 39 24 


70 


70 


-12.04 


3.66 


7.4 C 


IC2574 


UGC5666,DD081 


10 28 22.4 


68 24 58 


432 


213 


-17.17 


3.81 


7.85 ± 0.14 a 


UGC5692 


DD082 


10 30 35.0 


70 37 10 


153 


106 


-14.44 


3.80 


7.95 ± 0.20 b 


KDG73 




10 52 55.3 


69 32 45 


63 


50 


-10.75 


4.03 




NGC3741 


UGC6572 


11 36 06.4 


45 17 07 


100 


54 


-13.01 


3.24 


7.62 ± 0.20 b 


NGC4163 


NGC4167,UGC7199 


12 12 08.9 


36 10 10 


106 


82 


-13.76 


2.87 


7.91 ± 0.20 b 


UGCA276 


DDO113,KDG90 


12 14 57.9 


36 13 08 


80 


72 


-11.61 


2.95 




UGCA292 




12 38 40.0 


32 46 00 


58 


38 


-11.36 


3.62 


7.30 ± 0.03 b 


UGC8091 


GR8 


12 58 40.4 


14 13 03 


62 


46 


-12.00 


2.08 


7.65 ± 0.06 b 


UGC8201 


DD0165 


13 06 26.8 


67 42 15 


131 


75 


-15.09 


4.57 


7.63 ± 0.08 a 7.80 ± 0.20 b 


UGC8508 




13 30 44.4 


54 54 36 


80 


60 


-12.95 


2.58 


7.89 ± 0.20 b 


UGC8651 


DD0181 


13 39 53.8 


40 44 21 


97 


70 


-12.94 


3.14 


7.85 ± 0.04 b 


UGC8760 


DD0183 


13 50 51.1 


38 01 16 


108 


56 


-13.08 


3.22 


7.6 C 


UGC8833 




13 54 48.7 


35 50 15 


60 


58 


-12.31 


3.08 




KKR03 




14 07 10.7 


35 03 37 


34 


21 


-8.49 


1.97 




UGC9128 


DD0187 


14 15 56.5 


23 03 19 


64 


11 


-12.43 


2.21 


7.75 ± 0.05 b 


UGC9240 


DDO190 


14 24 43.5 


44 31 33 


111 


91 


-14.14 


2.79 


7.95 ± 0.03 b 


KKH98 




23 45 34.0 


38 43 04 


63 


10 


-10.29 


2.54 





Note. — Properties of the Sample Dwarf galaxies— Columns 1 and 2: galaxy names;Column 3 and 4: J2000 Right Ascension 
(RA) and Declination (DEC) as reported by the NASA/IPAC Extragalactic Database (NED);Columns 5 and 6: semi-major and 
semi- minor axis measured by the 3.6 /an images by Dale et al. (2009);Column 7: absolute B magnitude ( Karachentscv ct al. 
2004);Column 8: adopted distance (Tip of the Red Giant Branch method) (Dalcanton et al. 2009);Column 9: metallicity taken 
from different sources: IMoustakas et al.l I pOlOT ). ' iMarble et"aT] l(2010h . ILee et al.l 1(20071 1. 



To use the ground-based and HST images as a uniform 
data set the HST images needed to be smoothed and re- 
binned to the ground-based images' resolution. This step 
is necessary to achieve consistent magnitudes and colors 
because the HST images contain less light from nearby 
contaminating sources. Note that all of the clusters are 
unresolved in our ground-based images. 

The smoothing was performed using the IDL rou- 
tine FILTERJMAGE which convolves the image using a 
Gaussian kernel. The FWHM of the Gaussian kernel was 
set to the typical seeing of the ground-based images of 
1//5. The re-binning was carried out via the IDL routine 
HASTROM and has the effect of enlarging the pixel size 
of space-based images to facilitate a comparison with the 
ground-based data. 

The cluster photometry was carried out in IDL with 
the routine APER, which was adapted from DAOPHOT. 
Each cluster's photometry parameters were set via visual 
inspection of radial plots from the IRAF task PHOT in 
the interactive mode. Clusters in non-crowded regions 



were given an aperture of twice the PSF FWHM and 
clusters in crowded regions were given aperture radii just 
below the radius at which contaminating light becomes 
prevalent. To measure the local background, annuli were 
set at least 2 pixels from the edge of the photometric 
aperture. The pixel scale of the ground-based images and 
the convolved HST images is 0.5 arcseconds per pixel. 

Aperture photometry was performed on at least 10 un- 
saturated field stars at the two apertures and then at 5 
times the FWHM. The average and standard deviation 
of the flux ratios were used as the aperture correction 
and its uncertainty. 

HST zcropoints were obtained from the Space Tele- 
scope Science Institute website^ After the recovery 
of the ACS camera on July 4, 2006 the operating 
temperature changed and resulted in a revised set of 
photometric zeropoints. The appropriate zeropoints 
used were based on the observation date of each im- 

6 http://www.stsci.edu/hst/acs/analysis/zeropoints 



age. The HST data were then transformed into the 
Jo hnson-Cousins magn itude system via the prescription 
of iSirianni et al.l (12005ft. We corrected for Galactic ex- 
tinction using the ISchlegel et al.l (J1998I ) maps tabulated 
by the NASA/IPAC Extragalactic Database^ 




U-B 



Fig. 1. — The color-color distribution of the cluster sample with 
ages less than 100 Myr . The shaded areas indicate the density of 
the stochastic models covering ages less than 100 Myr with Z=0.004 
(SMC). The models lie primarily on a locus, but with a plume 
to red colors where red super giants have evolved off the main 
sequence. 



2.4. Cluster Properties 

Cluster ages and masses were determined usi ng the 
Bayesian method of iFouesneau fc Lancon (2010). The 
method is based on a large collection of Monte-Carlo 
simulations of individual clusters, each containing a 
finite number of stars, which explicitly accounts for 
the stochastic effects induced by the random presence 
of luminous stars. The synthetic clusters are con- 
structed with the population syn thesis code PEGASE.2 
(|Fioc fc Rocca-Volmerangel 119991) . The underlying stel- 
lar evolution tracks are those of the Padova group 
(jBressan et alJll993l ). with a simple extension through 
thermally pulsating Asymptot ic Giant Branch (AGB) 
based on the prescriptions of iGroenewegen fc de Jona 
(1993). The input stellar spe ctra are taken from the li- 
brary of iLeieune et al.1 (|1997l ) . The stellar I nitial Mass 
Function (IMF) is taken from iKroupa et ail (|1993t ). and 
extends from 0.1 to 120 M Q . Nebular emission (lines 
and continuum) is included in the spectra and broadband 
fluxes under the assumption that no ionizing photons es- 
cape. 

The method i s ex plained in detail by 
IFouesneau fc Lancon! (J20 1 0) , we therefore only briefly 
summarize the process here. Discrete model clusters 
will have a range of luminosities for a given mass due 
to the stochastic population of luminous stars. Hence 
mass-to-light ratios can not be treated as a constant for 
a cluster of a given age. For this reason, the total mass 

7 http://ned.ipac.caltech.edu/ 



of a cluster cannot be determined by a simple scaling 
factor of the total light, and must be simultaneously 
fit along with other cluster parameters. The Bayesian 
method establishes joint probability distributions of 
the intrinsic cluster parameters (age, mass, metallicity, 
extinction), given a set of photometric measurements 
and their uncertainties. Each intrinsic property of a 
cluster is estimated from the comparison of its whole 
set of absolute fluxes to the distribution of a model 
catalog. The set of intrinsic properties with the largest 
probability is then assigned to the cluster. 

Although this method is more advanced than a \ 2 min- 
imization, it still incorporates multiple free parameters 
(e.g., age, mass, metallicity, internal extinction, etc). In- 
dependently quantifying some of these parameters or us- 
ing reasonable assumptions can greatly increase the ac- 
curacy of derived cluster properties. We make two sim- 
plifying assumptions: (1) that the extinction is negligi- 
ble in these low mass dwarf galaxies and (2) we fix the 
metallicity based on estimates of the current gas phase 
metallicity. 

The lack of internal extinction in dwarfs has been 
demonstrated by many local volume surveys, which 
find lower infrared t o ultraviolet rati o s for a given ul- 
traviolet color (e.g., [Buat et al.l [2005t iDale et all 120071: 






iMunoz-Mateos et al.|[2009t: iDale et all |2009j). Addition 
ally, the well defined main sequence and blue helium 
burning populations seen in the CMDs of these galax- 
ies show that the interna l reddening is small ev en for 
young stellar populations ()Weisz et alJl2008l l201lT ). 

Fortunately, the metallicity of many of our galaxies 
have been observationally determined via spectroscopic 
studi e s of HII regions; listed in Table 1 (jMarble et al.l 
120101 : iMoustakas et al.1 120101) The current available 
models of IFouesneau fc Lanconl (|2010l ) have SMC, LMC, 
and Solar values (log(O/H)+12~8.0, 8.3, and 8.69, re- 
spectively) available at this time. The majority of our 
galaxies with clusters show similar metallicities to that of 
the SMC, however two galaxies (KDG61 and UGC5336) 
have values similar to the LMC. There are five galax- 
ies (DDO 78, IKN, UGC5442, UGC5428, & UGC8760) 
with clusters that do not have reliable metallicity mea- 
surements: these are fit with an SMC model. Using 
the Luminosity-Mctallicity relationship of nearby galax- 
ies (z < 0.1) we have estimated that all five of these 
galaxies have a metallicity of 12-r-log( Q/H)^7.5, given 
the a bsolute B magnitude of the galaxy ()Lee et al.ll2003L 
I2007D . The mass in clusters of these galaxies will be bi- 
ased to somewhat higher masses by approximately 5% 
given the Padova SSP models. 

The age ranges we choose to study are 4-10 Myr 
(<10 Myr) and 4-100 Myr (<100 Myr). We selected 
these age ranges to facilitate a direct comparison be- 
tween our data and those of previous studies of cluster- 
host relationships. Furthermore, these age ranges are 
characteristic of the time scales of the Ha and FUV 
SFR indicators (<10 and <100 Myr, respectively). The 
sh orter age range (<10 Myr ) mimics the age ranges 
of iGoddard et al.1 (|2010r) and (JAdamo et al.l 1201 ID while 
the longer age range (<10 Myr) mimics those used by 
iSilva- Villa fc Larsenl (|2011l) . The Padova isochrones used 
to determine the SFHs are not computed below an age 
of 4 Myr and are therefore not included in the SFHs of 
iWeisz et al.l (|2011l) . However, the isochrones used to de- 



termine the cluster parameters do extend below 4 Myr 
and our formal cluster age dating extends to clusters of 
1 Myr. Despite the differences in formal age limits, we 
expect the SFHs and the cluster ages to probe the same 
time bins. 

We determine the approximate completeness limit of 
our sample by examining the V magnitude histogram 
of all clusters. The peak in the number of clusters at 
the magnitude bin of 21.5 is used as the completeness 
limit. Our completeness limiting magnitude corresponds 
to ^1500 and ^3000 M Q at the median distance of our 
sample for the two age ranges we study in this analysis 
of <10 and <100 Myr age ranges, respectively. Although 
the peak of clusters in the V magnitude histogram is a 
rough indication of our completeness limit, we will show 
in Section 5.5 that the limit does not significantly change 
the overall results of our study. 

Figure 1 shows the colors of the clusters and the 
stochastic models over-plotted as a shaded region, where 
the ages of clusters and models are less then 100 Myr. 
Figure 2 shows the inferred age-mass diagram for all the 
clusters with an age less than 1 Gyr. The clusters with 
ages to the left of the vertical lines at log (age) of 7 and 
8 constitute the cluster sample of the two age ranges 
(<10 Myr and <100 Myr) used in this study. It is inter- 
esting to note the lack of massive (M > 4.5 x 10 4 M©) 
clusters in our selected age ranges compared to older 
ages. 
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Fig. 2. — The age-mass diagram for all the clusters with an age 
less than 1 Gyr. The vertical lines at log (age) of 7 and 8 repre- 
sent the upper end of the two age ranges used in this study. The 
horizontal arrows are visual representations of the two age ranges, 
<10 and <100 Myr. 



sample. 

We fit the data above the completeness limit and find 
good agreement with multiple studie s that show a nearly 
uniform power-law slope of —2 ± 0.2 (Zhang & Fall 1999; 
i Hunter et afl[2 003: Bi k et all 120031 Ide Griis et all 2003: 
McCradv fc Graham! I2007t iFouesneau et alJl2012f >. The 
X 2 fit yields a slope of 1.94±0.26. We also fit a power-law 
to the data with a fixed slope of -2 to illustrate the dif- 
ference between the fitted slope and the canonical slope. 
These slopes are graphically displayed in Figure 3 as solid 
and dashed lines for the fitted and fixed slopes, respec- 
tively. More recent studies show evidence that the CMF 
is more accurate ly represented as a Schechter function 
(|Schechterlll976l) . where the low cluster mass regime is 
characterized as a power-law and the high mass regime 
exponentially decays (M* > 10 5 Mm; iGieles et al. 2006 
Bastianl [2001 iGieied wm lLarseiJ 1200a iBastian et"al 



20121) . However, our single component power-law fit is 



justified since none of our clusters are near this high mass 
elbow in the Schechter function. 



fixed slope = -2.0 

fitted slope =-1.94±0.26 

<100 Myrs Clusters 



-2.0 - 



\. -2.5 - 



-3.0 - 




log (M/M a ) 

Fig. 3. — The Cluster Mass Function (CMF) of the cluster sample 
in the <100 Myr age range. The vertical error bars represent the 
Poisson error and the horizontal error bars show the mass range 
covered by each mass bin. The vertical dotted line represents the 
completeness mass at the median distance of the sample computed 
from the discrete models. Starting at the completeness limit, we 
populate equal number of cluster bins (9 per bin) to avoid small 
number biases. We also populate one bin below the completeness 
limit with all clusters below the completeness limit; this bin is 
plotted as an open square. The solid line is the fit to the data 
above the completeness limit. The dashed line is the fit with a fixed 
slope of -2 to facilitate a visual comparison between the canonical 
CMF slope and the unconstrained fit. The fit to the data shows 
agreement with a canonical CMF slope of -2. 



Figure 3 shows the Cluster Mass Function (CMF) for 
the clusters with an age less then 100 Myr. We populate 
mass bins with an equal number of clusters (nine per bin) 

to av oid small number biases (|Mafz Apcllani z~fe Ubedal 
120051 ) . We also populate one bin with all of the clusters 
below the completeness limit; this bin is plotted as an 
open square and shows a dearth of clusters. The verti- 
cal error bars represent the Poisson statistics of the bin 
and the horizontal error bars show the mass range for 
each bin. The vertical dotted line corresponds to our 
completeness limit at the median distance of the galaxy 



2.5. Star Formation Histories 

The Star Formation Histor y (SFH) of each galaxy 
came from iWeisz et all (1201 If ) and was g e nerate d us- 
ing MATCH (Dolphin 2002). Weisz et all (|201ll) used 
a Salpeter IMF (|Salpet cr 1955) with the Pad ova stel- 
lar evolutionary models of Mari go et all (l2008f ). These 
models are similar to those used for cluster fitting, 
but combines updat ed AGB tracks wi t h the models of 
iBertelli et ail (fl99l and iGirardi eFall (|2002T) . 

A full description of the method used to generate SFHs 
can be found in iDolphin] (J2002), however, we summa- 



() 



rize the key points here. The SFH code allows the user 
to input an IMF, binary fraction, and allowable ranges 
in age, metallicity, distance, and extinction. Synthetic 
Color Magnitude Diagrams (CMDs) are created for ev- 
ery combination of the input parameters and compared 
to the observed CMD using a parametrized % 2 . Mini- 
mizing the x 2 yields the most likely SFH of the galaxy. 
Due to the presence of luminous main sequence and su- 
pergiant stars, the SFHs over the l ast 100 Myr consid- 
ered can be derived very accurately (|Dohm-Palmer et al.l 
H99l IWeiTz et al.ll200ll . 

3. CLUSTER AND STAR FORMATION PROPERTIES 

In this section we describe the calculation the cluster 
formation efficiency (r — CFR/SFR) and star forma- 
tion rate density (T*sfr)- We then present the raw data 

and examine the T-Y,sfr and M v Ilg 1 os -SFR relations 
in our dataset. As we will show, no clear trend is seen 
in the T-'Sgpfi relation, while the My i& cs -SFR shows 
the expected correlation. However, both plots show sig- 
nificant scatter whose source will be discussed in depth 
in Section 5. 

3.1. Cluster Formation Rates 

The cluster formation rate of each galaxy is computed 
by adding up the mass of clusters in each age range then 
dividing by the duration; these values are listed in Table 
2. We have made a power- law CMF correction for each 
galaxy to account for missing faint clusters. The slope 
of the power- law is set to —2 and we have made the 
assumption that the least massive bound cluster would 
have M > 100M Q . The assumption of a pure power- 
law CMF is valid due to the lack of any massive clusters 
and the agreement of the sample's CMF with the canon- 
ical —2 power-law (see Figure 3). The missing mass for 
galaxies with cluster detections are computed by inte- 
grating under the CMF from 100M Q to the lowest mass 
bin where clusters are detected. Since we verified that 
each individual galaxy's CMF showed no significant de- 
viation from a canonical CMF, we did not integrate to 
the completeness limit of each galaxy. 

There are some dwarf galaxies in our sample that do 
not have any detectable clusters in the age ranges stud- 
ied. These galaxies will be represented as an upper limit 
symbol (downward arrow). In the Y-Y,sfr figures the up- 
per limits correspond to the maximum total cluster mass 
undetectable and are calculated by integrating under a 
CMF with a slope of —2 between 100 Mq and the lim- 
iting completeness mass for each galaxy. The complete- 
ness mass is determined from our completeness magni- 
tude limit (my =21.5). At the distance of each galaxy 
the completeness limiting magnitude is translated into a 
mass given the Bayesian model median mass for each age 
range. Since there are no clusters detected just above the 
completeness limit in galaxies without detected clusters, 
we assume that at most one cluster will be present at the 
completeness limit. This assumption translates into the 
amplitude of the CMF for galaxies without detected clus- 
ters. Dividing the missing mass under the CMF by the 
duration of each age range yields an upper limit CFR. 
Dividing the limiting CFR by the measured SFR pro- 
duces an upper limit cluster formation efficiency (r). In 

the M v llg tcst -SFR figures the upper limits correspond 



to our completeness V magnitude of 21.5 at the distance 
of each galaxy. 

We determine the cluster formation rates in two dif- 
ferent age ranges, <10 Myr and <100 Myr. The cluster 
formation rates in the <10 Myr age range are based on 
small numbers of clusters and also suffers from uncer- 
tainties in the age. However, we have verified that each 
of these clusters show moderate to strong Ha emission 
which traces star formation of the appropriate age range 
(see Section 3.3). The <100 Myr age range clusters will 
have larger number statistics due to the longer duration, 
but here the effects of cluster destruction may be impor- 
tant. 

An IMF correction needs to be applied to our cluster 
sample due to the different adopted IMFs assumed for 
cluster parameter (Kroupa IMF) and the SFR (Salpeter 
IMF) determinations. A Kroupa IMF deviates from a 
Salpeter IMF to lower numbers of stars at low masses 
and therefore underestimates the mass in clusters if the 
CFR and SFR are compared. Integrating under both 
IMFs normalized to produce the same number of stars 
yields a mass ratio of 1.38. We multiply our CFRs by 
this factor to account for the different IMF assumptions. 

3.2. Star Formation Rates 

We calculate the field SFR from the CMD-based SFHs 
(e. g., the SFR as a function of time and metallicity) 
of IWeisz et"aLl (|2011[ ). The SFHs were generated with 
the method described in Section 2.5. The correspond- 
ing SFRs are calculated by taking the total stellar mass 
formed during the age range and divided by the dura- 
tion. The resulting average SFRs are listed in Table 2. 
It should be noted that there may be periods of higher 
or lower SFRs during the adopted time interval. 

3.3. T,sfr 

The star formation rate surface density (Ssfr) is cal- 
culated by normalizing the SFR by the area of the galaxy 
being studied: typically the whole galaxy in our sam- 
ple. Most of our galaxies have a small enough angu- 
lar size to fit entirely within HST ACS camera field of 
view. The size of these galaxies were determined us- 
ing th e area defined by the 3.6 (im ellipses of Dale et al. 
(2009), listed in Table 1. Since FUV and Ha flux are SFR 
indicators of 100 Myr and 10 Myr star formation, respec- 
tively, we have visually verified that these ellipses con- 
tain both the FUV and Ha emission in each galaxy. The 
FUV and Ha images were taken as part of the 11 Mpc 
Ha UV Galaxy Survey CllHUGSi lKennicutt et al"]l200i 
ILee et all2009al) In five galaxies, the 3.6 /xm ellipses were 
larger than the ACS FOV. The normalizing area for these 
galaxies were taken as the number of pointings times the 
ACS FOV (202" x 202"). 

3.4. T-Zsfr 

With both the cluster formation rates and the star for- 
mation rates calculated, the cluster formation efficiency 
can be quantified. By definition V is the ratio of the CFR 
divided by the SFR. 

Figure 4 shows the T-T, S fr plot. The <10 Myr data 
show a large scatter to higher Ts, whereas the <100 Myr 
data show lower T values. The downward arrows are 
galaxies where no clusters were detected and represent 



TABLE 2 
Dwarf Galaxy Clustered Star Formation 



Galaxy 




(4-10 Myr) 








(4-100 Myr) 








(CFR) 1Q 


CMF Corr 


(SFR) 10 


7i j-brightcst 

(mag) 


(CFR) m0 


CMF Corr 


(SFR) 100 


7i ^brightest 

(mag) 




(M /yr) 


(factor) 


(M /yr) 


(M Q /yr) 


(factor) 


(M Q /yr) 


ESO540-G030 






0.0 








2.28c-05 




NGC404 


















KKH37 






0.0 








4.71e-04 




NGC2366 


4.07e-04 


3.29 


9.35e-02 


-8.52±0.02 


1.73e-04 


1.86 


6.26e-02 


-8.52±0.02 


UGCA133 






0.0 








1.87e-05 




UGC4305 


3.01e-03 


5.04 


1.14e-01 


-8.88±0.01 


7.51e-04 


2.26 


6.37e-02 


-8.88±0.01 


M81da 






2.07e-03 








1.66e-03 




UGC4459 


5.60e-04 


7.94 


3.25e-03 


-7.90±0.03 


4.84e-05 


5.75 


6.53e-03 


-7.90±0.03 


UGC5139 






1.75e-02 








1.43e-02 




FM2000.1 






3.59e-06 








6.00e-06 




BK3N 






1.36e-05 








1.57e-03 




KDG61 






4.71e-04 








1.02e-04 




arpsloop 






0.0 








2.54e-03 




UGC5336 


9.86e-04 


2.76 


3.71e-03 


-8.59±0.01 


2.54e-04 


2.01 


2.69e-02 


-8.59±0.01 


LEDA166101 






0.0 








1.34e-05 




UGC5428 






0.0 








0.0 




UGC5442 


















IKN 






0.0 








1.64e-04 




HS98.117 






1.42e-04 








8.91e-06 




DD078 






5.20e-06 








2.18e-05 




IC2574 


4.10e-03 


2.10 


1.06e-01 


-9.12±0.03 


1.29e-03 


1.55 


8.18e-02 


-9.12±0.03 


UGC5692 


2.14e-04 


5.00 


7.34e-03 


-9.44±0.03 


1.85e-05 


3.62 


1.28e-03 


-9.44±0.03 


KDG73 






0.0 








8.14e-04 




NGC3741 






1.44e-02 








4.22e-03 




NGC4163 






9.93e-03 








2.58e-03 




UGCA276 






0.0 








1.98e-05 




UGCA292 






5.15e-03 




1.26e-04 


2.82 


3.35e-03 


-7.46±0.04 


UGC8091 






1.65e-02 








1.77e-03 




UGC8201 






3.16e-02 




1.29e-03 


1.32 


5.79e-02 


-9.16±0.02 


UGC8508 






6.54e-03 








2.77e-03 




UGC8651 






8.34e-04 








3.41e-03 




UGC8760 


2.14e-04 


5.00 


3.03e-03 


-7.37±0.08 


1.85e-05 


3.62 


4.29e-03 


-7.37±0.08 


UGC8833 






2.52e-03 








1.56e-03 




KKR03 






0.0 








2.42e-04 




UGC9128 






2.79e-04 




1.77e-05 


3.74 


1.29e-03 


-5.45±0.12 


UGC9240 






5.24e-03 




6.68e-05 


4.44 


8.86e-03 


-5.90±0.19 


KKH98 






0.0 








4.97e-04 





Note. — The cluster formation and star formation results organized by two age ranges (<10 Myr and <100 Myr), 
as indicated by the top label. For each age range this table presents the Cluster Formation Rate (CFR), the Cluster 
Mass Function (CMF) correction used in Section 4, the CMD-based Star Formation Rate (SFR), and the V absolute 

magnitude of the brightest cluster (M„ ng os ). The CFRs include the IMF correction factor of 1.38 presented in Section 



3.1. An ellipsis represents no data or no detectable clusters. 

upper limits on the cluster formation efficiency for these 
galaxies. Low mass clusters below our completeness limit 
may be present, but the total cluster mass will fall below 
the upper limit arrows. 

Due to the small numbers of clusters in each of our 
dwarf galaxies, we expect the scatter seen in our data to 
be dominated by stochastic effects. The propagation of 
formal errors in cluster age and mass does not include 
this scatter. In Section 5 we model the stochastic effects 
and provide resulting confidence intervals for each galaxy. 



3.5. M. 



brightest 
V 



SFR 



The scatter in T is motivation for analyzing another pa- 
rameter space that compares cluster formation and star 
formation, namely the V — band luminosity of the bright- 
est cluster (M v " s cs ) versus the star formation rate. 
The brightest cluster in a galaxy has previously been 
found to correlate with the galaxy's SFR, as a larger SFR 
results in a la rger cluster p opulation sampling higher 
mass clusters (Larsen 2002). In Figure 5 we plot the 
brightest cluster in the appropriate age range (<10 Myr 



or <100 Myr) versus the galaxy's SFR. Both age ranges 
show a positive trend with a moderate amount of scat- 
ter. The upper limit data points are galaxies with no de- 
tectable clusters and represent the absolute V magnitude 
of our completeness limit at the distance of each galaxy. 
At the higher SFR end of our sample, these upper limits 
are often 3-4 magnitudes fainter than galaxies with de- 
tected brightest clusters and similar SFRs. Stochasticity 
in the number of rare, massive clusters may contribute 
significantly to the scatter; we model the stochastic ef- 
fects on this relationship in Section 5. 

4. RESULTS 

In this section we discuss the results of previous stud- 
ies on the relation between star formation and cluster 
formation in higher mass galaxies. We take these stud- 
ies' results as presented in the original work, highlighting 
sources of uncertainty and differences from out analy- 
sis. We then present a comparison of our T-Y>sfr and 
My lg ' cs -SFR relations to previous work done in higher 
mass galaxies. In Section 5 we present a detailed analysis 
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Fig. 4. — The cluster formation efficiency (r=CFR/SFR) versus 
the star formation rate density (Eg^^^SFR/kpc 2 ) for two age 
ranges. The data have been corrected for incompleteness via a 
CMF with a slope of -2. Data from the <10 Myr interval are 
presented as red asterisks and the <100 interval as blue X's. The 
downward arrows are galaxies where no clusters were detected and 
represent upper limits on the total mass in clusters; the upper limits 
have the same color coding as the non-upper limit data points. 
The <100 Myr age data show up to a factor of 100 variation in the 
cluster formation efficiency at fixed star formation rate density. A 
large amount of this scatter is toward low cluster formation rates. 
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Fig. 5.— The M^ rightost -SFR for the two age ranges. Data from 
the <10 Myr interval are presented as red asterisks and the <100 
interval as the blue X's. The downward arrows are galaxies where 
no clusters were detected and represent upper limits on the total 
mass in clusters; the upper limits have the same color coding as the 
non-upper limit data points. Both age range's brightest clusters 
show a positive trend with SFR. 



of the scatter in these relationships using simulations of 
stochastic cluster formation. 



TABLE 3 
Other Studies Results 



Galaxy 



CFR 
(M Q /yr) 



SFR 
(M /yr) 



Area 
(kpc 2 ) 



r 

(%) 



Rcf. 



NGC1569 


0.05 


0.3626 


13.0 


13.9±0.8 


a 


NGC3256 


10.57 


46.17 


74.9 


22 Q+ 7 - 3 
^ ,a -9.8 


a 


NGC5236 


0.103 


0.3867 


0.708 


26 7+ 5 - 3 


a 


NGC6946 


0.022 


0.1725 


37.5 


12.5±I;1 


a 


LMC 


0.007 


0.1201 


79.0 


5.8±0.5 


a 


SMC 


0.002 


0.0426 


58.9 


4 2+ ' 2 


a,d 


Milky Way 


0.011 


0.1508 


12.6 


7 0+ 7 - 


a 


ESO338-IG04 


1.6 


3.2 


2.07 


50±10 


b 


Haro 11 


11.2 


22.0 


10.2 


50±g 


b 


ES0185-IG13 


1.7 


6.4 


12.3 


26±5 


b 


Mrk 930 


1.33 


5.34 


9.05 


25±10 


b 


SBS0335-052E 


0.64 


1.3 


1.37 


49±15 


b 


NGC5236 


0.038 


0.39 


28.7 


9.8 


c 


NGC7793 


0.015 


0.15 


23.1 


9.8 


c 


NGC1313 


0.061 


0.68 


60.0 


9.0 


c 


NGC0045 


0.009 


0.05 


49.0 


17.3 


c 


NGC4395 


0.005 


0.17 


36.5 


2.6 


c 



Note. — The results from previous studies presented here 
are from: " IGoddard eTal! II20RT) . ^Adamo et~ al. (2011), and 
ISilva- Villa fc Lar son (2011), and ^iclcs & Bastian (2008]). The 
quantities listed are the galaxy name, Cluster Formation Rate 
(CFR), the Star Fo rmation Rate (SFR), the normalizing are a, and 
T (CFR/SFR). The lGoddard et al.l l[2010l ') and lAdamo et alJ pOTl 
values have the IMF correction factor of 1.38 incorporated into 
them. 



end. Both iGoddard et al.1 (|2010D and lAdamo et all 
(120111) used SFRs der i ved fr om integrated fluxes, whereas 
ISilva- Villa fe Larsenl (I20TI used CMD-based SFRs de- 
rived from a similar method as the one presented in 
this study. The integrated flux SFRs were based on 
the K ennicutt-Schmidt law which assumes a Salpeter 
IMF (Kcnnicutt 1998). To achieve consistency be- 
tween the IMFs used in estimating the cl uster formation 
rate, t he clu ster masses from this study, IGoddard et al.1 
(J2010D . and lAdamo et all (120111) needed to be multi - 
plied by a factor of 1.38. ISilva- Villa fc Larsenl (|2011[ ) 
used the same IMF in determining both quantities. The 
characteri stic time scales for t hese studies a r e <10 Myr 
for IGoddard et al.l (I20T0I) andlAdamo et al.l (|2011[ ). and 
<100 Myr for ISilva- Villa fc Larsenl (120 111). 

Add itional efforts were made by ISilva- Villa fc Larsenl 
(|2011l ) to investigate the effects of cluster destruction 
on r by modeling the destruction-dependent luminos- 
ity function. However, we treat the T versus T,$fr 
relationship as purely observational and therefore only 
use their values labeled P I (Paper I analysis method; 
ISilva- Villa fc Larsenl f2010h . The PI values are those ob- 
tained without cluster destruction taken into account. 
The results of all three studies are listed in Table 3. 



4.1. Previous Studies 

All of the previous studies of cluster formation effi- 
ciency share several common elements. In each, clus- 
ter masses were determined using continuous models, 
as opposed to the discrete models used in our anal- 
ysis. Also, each uses a cluster mass function correc- 
tion assuming a —2 slope power-law at the low mass 



4.2 






We combine our T-T<sf r me asurements wit h 
data from the studies of IGoddar d et al.l ( 2010), 
ISilva- Villa fc Larsenl ()201lD . and lAdamo et al.1 (|2011l ) in 
Figure 6. The solid line is the linear regression fit to the 
data of IGoddard et all (|2010f) . The dotted lines above 
and below the trend line are three times the standard 



deviation of all the previous study's d ata around the 
trend line. Although, the relationship of Goddar d et al.l 
(2010) is supported by many sources, none of them are 
probing the star formation rate regime of our data. 
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Fig. 6. — The FSsfr plot, similar to Figure 4, but with the 
previous study's d ata over-plotted as op en symbols. The solid line 
is the trend from Goddard et al. (2010) and the dotted lines are 
the 3 times the standard deviations of all the previous data. Data 
from both age bins show poor agreement with the trend of previous 
studies. The majority of the galaxies in the <100 Myr age range 
fall significantly below the trend indicating a dearth of total cluster 
mass. To study the effects of low number statistics we have binned 
our data for the two age ranges. The <10 Myr and <100 Myr 
binned data points are the large, red, open star and large, blue, 
open X, respectively; the x-axis error bars represent the width of 
the bin. The binned <10 Myr data show good agreement with the 
expected trend, but the <100 Myr data show a dearth of clusters 
despite the better number statistics of the longer age range. 



We find poor agre ement for both age ra nges with the 
extrapolation of the iGoddard et al.l (|2010l) relationship. 
The <10 Myr age range data are discrepant to both high 
and low T values whereas the majority of the <100 Myr 
age range galaxies show a dearth of cluster formation. 
UGC5692 deviates above the trend line in both age 
ranges and NGC2366 deviates below the line in both age 
ranges. This scatter suggests there may be a true varia- 
tion in the cluster formation properties of these galaxies. 
We discuss the outliers in more detail in Section 5. 

To examine the effects of low number statistics we have 
binned all of the clusters into one cluster population. We 
performed an analysis of the combined cluster popula- 
tion by treating these clusters as if they occupy a single 
galaxy. The procedure is the same as described in Sec- 
tion 3.1. We calculate the total cluster mass in the <10 
Myr and <100 Myr age ranges, including all galaxies with 
\og(T,sFR) values > —4.5. The binned T values are cal- 
culated by dividing a CMF corrected cluster formation 
rate by the sum of the star formation rates; including the 
upper limit star formation rates. The ^sfr values are 
determined by summing the SFRs and dividing by the 
sum of the areas. The binned <10 and <100 Myr data 
are shown on Figure 6 as a large, red, open star and a 
large, blue, open X, respectively. 

The <10 Myr data show good agreem ent with the ex- 
pected trend of IGoddard et al.l (|2010l ). However, the 
binned <100 Myr data shows a significant deviation to 



lower r values (lower total cluster mass). Due to its 
longer age-range, the <100 Myr age range has more clus- 
ters and better number statistics than the <10 Myr data. 
However, it contains fewer clusters than expected from 
the T-T,sfr relationship. We will discuss possible sources 
of the dearth of clusters in the <100 Myr age range in 
Section 5. 



4.3. 



M^ ightcst -SFR 



Plotting the brightest cluster versus the galaxy-wide 
SFR facilitates the analysis of clust ered star form ation 
by directly testing stochastic effects. iLarsenl (|2002t ) orig- 
inally proposed the statistical "size-of-sample" effect ex- 
planation for this correlation, such that higher SFR 
events can sample higher into the CMF to create higher 
mass (assumed brighter) clusters. In other words, this 
relationship holds because of stochastic effects. Figure 7 
shows t he data from this study over-plott ed onto those o f 
ILarsenl (pOOl ). iBastianl (|2008| ). and I Adamo et al.1 (|2011f l. 
The trend lin e is a linear regre ssion fit of the dat a from 
ILarsenl (J2002H performed bv lWeidner et al.l (|2004t) . The 
dotted lines above and below the trend represent three 
times the standard deviation of all previous data points. 
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similar to Figure 5, but with the previous data over-plotted as 
o pen sym b ols. T he solid line is the linear regression fit to the data 
of lLarsenl J2002T) performed bv lWeidner et ail 1 120041) . The dotted 
lines are the 3 times the standard deviations of all the previous 
data. Both age ranges show good agreement with the established 
trend, but there are a few outliers. Upper limits (no detectable 
clusters) below the trend line are significant since these galaxies 
should have formed a detectable bright cluster, given the star for- 
mation rate. 



There are outliers falling both above and below the 
expected M^ rightost -SFR relation. The typical V mag- 
nitude error for data presented in Figure 7 is approxi- 
mately 0.03 magnitudes. The five upper limit low outliers 
near a log SFR(M G yr -1 ) of -2 belong to the <10 Myr 
age range. These upper limits are significant because 
the brightest clusters predicted by the trend are brighter 
than the completeness limits of these galaxies, but no 
clusters of these magnitudes are detected. Conversely, 
the <100 Myr data show no significant outliers below 
the trend. It is interesting to note that all of the T-T,sfr 
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low outliers show good to moderate agreement with the 

M brightcst_ SFR trend 

There is also one high outlier galaxy (UGC5692) in 
the <100 Myr age range that lies significantly above the 
trend line. This galaxy is a lso a high outl ier in the T- 
^sfr relationship. However . IBastianl (|2008l ) showed that 
outliers above the trend were in fact not outliers after 
using the youngest (t agc < 10 Myr) clusters to show 
that the brightest cluster is more likely to be young 
and therefore more accurately reflects the recent SFR. 
This is true for UGC5692, whose brightest cluster has 
an age less than 10 Myr. UGC5692 has undergone a re- 
cent increase in SFR; using the <10 Myr data point (log 
SFR(Mq yr _1 )~ —3 and the same magnitude) brings 
this galaxy within the scatter of the relationship. How- 
ever, the lower limits at high SFRs in the <10 Myr sam- 
ple do see m to vio l ate th e M v " 6 cs -SFR correlation pro- 
posed by IBastianl (120081 ) . A full discussion of these out- 
liers be presented below. 

5. DISCUSSION: IS THE SCATTER DUE TO 
STOCHASTICITY 

The previous section showed that clustered star for- 
mation in dwarf galaxies does not necessarily agree with 
established cluster-host relationships. The sample shows 
a significant amount of scatter around these relationships 
with a large number of significant outliers, particularly 
in the T-Hsfr relationship. However, because the over- 
all SFRs of our dwarf galaxy sample are lower (due to 
both smaller galaxy mass and overall star formation ef- 
ficiency), the cluster population is much more likely to 
be affected by stochastic sampling of the cluster mass 
function. Therefore we examine the effect of stochastic 
sampling of the CMF, and whether this sampling can ac- 
count for the observed scatter in both the T-T,sfr and 

M brighte S t_ SFR relations _ 

5.1. SLUG 

We utilize the publicly availab le code "SLUG" (Sta- 
tistically Lighting Up Galaxies; Fumagalli e t al.l 120111 : 
Ida Silva et alJ feOll) to model these effects. SLUG cre- 
ates a simulated galaxy by generating clusters and stars 
via stochastic sampling of both the CMF and the stellar 
IMF, respectively. The code allows the user to turn on 
and off cluster destruction and to set the fraction of stars 
that form in clusters. 

In addition to providing a tool to incorporate the vari- 
ability in our data due to sampling of small numbers 
of clusters, the SLUG simulations also allow us to di- 
rectly compare our simulations to data without making 
a CMF correction for clusters fainter than our complete- 
ness limit. Half of the 17 galaxies with detected clusters 
have fewer than three clusters. Consequently, the CMF 
correction for these galaxies is highly uncertain as a re- 
sult of fitting the amplitude to one or two data points. 
The direct comparison enabled by SLUG eliminates this 
correction, and incorporates the uncertainty in this mass 
correction into our stochasticity estimates. 

5.2. t-z S fr 

To evaluate the scatter in T we use SLUG to simulate 
only the cluster population of each galaxy. The low val- 
ues of r of our galaxies mean that the amount of mass 



in clusters is a small fraction of the total star formation 
in the galaxies. Therefore, we expect the stochasticity 
of the cluster formation process to dominate the scatter 
in r, and simulate only the effects of stochasticity in the 
cluster population. 

For our simulations, we assume that a universal clus- 
ter formation mechan ism is accurately described by the 
iGoddard et~"aTI (|2010D relationship. We input into SLUG 
the CFR expected if each galaxy fell on the T versus 
^SFR relationship found bv lGoddard et al.l (|2010( 1. given 
the galaxy's measured T<sfr- SLUG allows the user to 
set the fraction of stars that form in clusters via a param- 
eter with values that range from < / < 1. Simulations 
can generate stars without clusters (/=0), clusters with- 
out stars (/=1), or some combination of both. Since we 
only want SLUG to simulate the cluster population of a 
galaxy, we input a CFR instead of the galaxy's SFR, we 
force all stars to form in clusters, and we disable cluster 
destruction. The combination of these input parameters 
effectively creates a cluster population based on the mass 
in clusters expected from the cluster formation rate ob- 
served in each galaxy. 

For each galaxy, we run 1000 simulations where each 
iteration randomly generates a new cluster population 
(ignoring field stars) from a CMF with a slope of —2 and 
a Salpeter IMF is assumed to populate each cluster. Our 
simulations also include the same stellar isochrones as 
those used for cluster fitting. The mass of clusters gen- 
erated with a luminosity brighter than our completeness 
limit is added up and divided by the measured SFR to 
produce a simulated I\ The median of the simulations 
is used as the prediction of the T — Y^sfr simulations, 
while the middle 68.3% of the simulations are taken as 
the 1<7 confidence interval. To make a direct comparison 
between the SLUG simulations and our observed cluster 
sample, we have also cut any real clusters fainter than 
our completeness limit when comparing our data to the 
simulations. 

While the absolute T values o f the simulations are 
based on the lGoddard et all ()2010f ) relationship, the scat- 
ter seen in the data points should generically reflect 
how much of the scatter in our data results from clus- 
ter stochasticity regardless of the relationship between T 
and T, S fr- 

The results for the <10 Myr and <100 Myr age ranges 
are graphically shown in Figure 8, in the top and bot- 
tom panels, respectively. Each galaxy is evenly spaced 
on the £-axis, which has been sorted by TIsfr with in- 
creasing values to the right. The open triangles are the 
SLUG simulated galaxies, and the red asterisks and blue 
x's are the <10 Myr and <100 Myr age range data points 
with detected clusters above our completeness limit, re- 
spectively. The simulated T values are set to a value 
of when the median of the simulations produced no 
clusters. Simulations for the majority of the lowest SFR 
galaxies (left of UGC5139 for the both age ranges) pro- 
duced very few detectable clusters in agreement with the 
lack of clusters observed in those galaxies. 

There are three major outliers (UGC5336, UGC8760, 
and UGC4459) to high T values in the <10 Myr age 
range, and one (UGC5692) in the <100 Myr age range. 
These galaxies fall noticeably above the simulations. 
None of the simulations for UGC5336 nor UGC4459 and 
less than 1% and 7% of the simulations for UGC8760 
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and UGC5692, respectively, produced T values as high 
as the observed data points. This makes UGC5692 a 
~ 2cr outlier and UGC8760 a - 3ct outlier. Since 1000 
simulations were run, the lack of overlap with UGC5336 
and UGC4459 make these galaxies > 3<r outliers. These 
galaxies provide strong evidence that more clusters are 
formed than expected from T-Hsfr relationship. 
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Also, there are large blue, nebulous active star form- 
ing regions in this galaxy that could hide young massive 
star clusters which may not yet be detectable. How- 
ever, we consider hidden clusters an unlikely explanation 
since many other galaxies in our sample (e.g., IC2574, 
UGC4305, and UGC4459) also contain large nebulous 
regions and all either have good agreement with or lie 
above the T-Y<sfr. trend in at least one age range. This 
implies that NGC2366 has truly formed fewer clusters 
than predicted by the T-~Esfr relationship. 

We have examined the clusters in the outlier galaxies 
carefully to insure that their discrepant T values are not 
due to erroneous age and mass determinations. Since our 
galaxies are close enough to partially resolve the stars in 
each cluster, we can broadly verify the age of the clus- 
ters in these galaxies. A young CMD will have a well 
defined main sequence and an old cluster will have a well 
defined red giant branch. All of the clusters in these 
outlier galaxies exhibit properties of a young stellar pop- 
ulation of stars within the them. Furthermore, we have 
double checked the discrepant <10 Myr clusters' ages by 
confirming moderate to strong Ha emission (SFR indica- 
tor of age less than 10 Myr) in each of these clusters and 
little to no Ha emission in older clusters. The high and 
low outliers strongly suggest that the cluster formation 
efficiency shows significant variation in dwarf galaxies. 

To test how well our data is described by the SLUG 
simulation model for the collection of galaxies as a whole, 
we calculate a reduced % 2 for each age range. We assign 
the observed data points as the expected value and treat 
the simulations as the "measured" data points with er- 
rors. We then evaluate the goodness of the fit by calcu- 
lating the reduced \ 2 between the simulations and the 
observed data: 



Increasing log Z SFR (M Q yr 'kpc 2 ) — > 

Fig. 8. — The simulated T results from SLUG, where the top 
and bottom panels are the <10 Myr and <100 Myr age ranges, 
respectively. Each galaxy is presented on the x-axis evenly spaced 
and sorted by Hgpji, with increasing values to the right. The 
open triangles are the SLUG simulated galaxies and the error bars 
are the 68.3% confidence inertvals of the fOOO realizations of the 
simulations. The blue X's are the observed data without CMF cor- 
rections and the downward arrows are the maximum total cluster 
mass for galaxies with no detectable clusters. The majority of the 
SLUG simulations in both age ranges agree with the observed data 
since they are within at least 2 sigma. However, there are a few 
significant outliers above and below the simulations which might 
indicate a true variation in the cluster formation of these galalxies. 
The <100 Myr age range simulations systematically fall above the 
observed data. This discrepency suggests that stochasticity cannot 
account for all of the scatter to low cluster formation seen in the 
<100 interval data. 



There is one galaxy that falls significantly below the 
T-Tjsfr trend line in both age ranges: NGG2366. This 
galaxy also falls nearly 3a below the SLUG simulatio ns in 
the <100 Myr age range as well. iBillett et ail (|2002[ ) also 
found this galaxy to have an unusual dearth of cluster 
formation compared to other galaxies with similar star 
formation properties. NGC2366 has spatially extended 
pockets of star formation where a lack of cluster forma- 
tion could be due to a lower gas density. The gas density 
of dwarf irregular galaxies at intermediate and extended 
radii is typically well below the thr e shold of l arge-scale 
star formation (Huntc r~et al.l [T998: Lcr ov et al.l [20081. 
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where N is the degrees of freedom (number of galaxies) , 
p is number of parameters being fit (one for the F val- 
ues), r sim is the median simulated data point, r orjs is the 
observed data point, and <7 s i m is the confidence interval 
of the simulations. Upper limits are included in the esti- 
mate of Xrod om y when they fall below the median of the 
simulations. The <10 Myr and <100 Myr x^ cd values 
are 3.6 and 7.2, respectively. 

Large \^ cd values could result from: (1) real deviation 
of the cluster formation efficiency from the T — T,sfr 
relationship, (2) cluster destruction; (3) non-Gaussian 
stochastic errors as determined from SLUG; or (4) ad- 
ditional errors in our data or analysis. Number four 
includes errors in the mass and age determination, the 
cluster identification, and the completeness limit. 

Most of the simulations in the <100 Myr age range 
show moderate agreement on a galaxy-by-galaxy basis 
with the observed data since they are within the 2ct 
spread of the SLUG simulations. However, the dwarf 
galaxy sample as a whole shows a systematic offset to 
lower r values compared to the stochastic simulations. 
The moderate agreement of most galaxies with the SLUG 
simulations in this age range imply that stochastic effects 
account for most of the scatter in T, but the systematic 
offset of the sample as a whole to lower T values and the 
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higher relative reduced % 2 value of the <100 Myr, com- 
pared to the <10 Myr age range, implies the presence of 
some additional source of scatter. 

One explanation for the lack of cluster formation in 
the <100 Myr age range, compared to the <10 Myr age 
range, is cluster destruction. The primary mechanism 
responsible for cluster destruction is currently under de- 
bate (see Section 5.3), but its effect should be significant 
in a longer age range (<100 Myr) compared to a shorter 
age range (<10 Myr). This effect is indeed observed in 
our data. Seven galaxies with cluster detections and four 
without detections (upper limits) in the <100 Myr age 
range show significant deviations from the T-T,sfr rela- 
tionship to lower r values (i.e. low total cluster mass), 
whereas only one galaxy with a cluster detection and one 
upper limit data point in the <10 Myr age range show 
a low r deviation (see Figure 6). In other words, 11 
data points in the <100 Myr age range are discrepant 
compared to only 2 in the <10 Myr age range. Further- 
more, the improved number statistics via binning the 
entire cluster sample shows that the binned <100 Myr 
data are also discrepent to lower T values compared to 
the binned <10 Myr data (see large, open symbols in Fig- 
ure 6). It is likely that some form of cluster destruction 
is responsible for at least some of the additional scatter 
in our sample. 

5.3. Constraints on Cluster Destruction 

More scatter towards low cluster mass is seen in the 
T — "Esfr relationship for the <100 Myr age range 
that cannot be fully accounted for by stochasticity. If 
cluster destruction is responsible for this deficit, then 
the My 16 cs -SFR relationship will provide evidence for 
whether or not the most massive clusters are surviv- 
ing. There are two popular models of cluster destruc- 
tion, (1) mass-independent destruction and (2) mass- 
dependent destruction. The mass-independent model 
is dominated by internal mechanisms, such as infant 
mortality, and does not preferentiall y destroy a clus- 
ter based on its' mass ilHillsl [l980t iFall etail 1 20051 
20091: iWhitmore et all [2007. 201fA iChandar et al .1 [2000. 
2010alrb1 ). The mass-dependent model, on the other 
hand, is dominated by external tidal forces such as 
proximity to giant molecular clouds or interior sec- 
tions of spiral galaxies, and thus preferentially destroys 
smaller mass clusters (Spitzcr & Harm 1958; Hcnon 



1961; Baumgardt & Makino 20 031: lLamers et al.l 120051: 
Gieles et al.l 120 07; Bast ian et al.l l2011h . Unfortunately, 
r does not distinguish which clusters are missing from 
the total mass of clustered star formation and there- 
fore cannot distinguish between these models. However, 
the My 1 ' 8 cs -SFR relation correlates the brightest (and 
therefore most massive) cluster to the SFR of the galaxy, 
and analysis of both relationships may provide evidence 
of a dominant cluster destruction model. 

The T-'Esfr relationship is sensitive to both cluster 
destruction models since both effect the overall cluster 
mass. However, the My" s cs -SFR relationship should 
be unaffected by a mass-dependent destruction model 
and should show no major deviation from the trend 
since a mass-dependent cluster destruction model pre- 
dicts that massive clusters will preferentially survive. 
If a mass-dependent model is dominant then the low 



r galaxies should show agreement with the M v " g es - 
SFR relationship. This cross-relationship prediction is 
seen in our data. Of the 11 galaxies with significantly 
low r values, 4 out of 4 upper limits and 7 out of 7 
galaxies with cluster detections show good to moder- 
ate agreement with M v Ilg es -SFR relation. This sug- 
gests that the missing mass in the cluster formation seen 
in the <100 Myr age range is not due to random de- 
struction, but rather mass-dependent cluster destruction 
where the high mass (brightest) clusters preferentially 
survive. This is evidence for mass-dependent cluster de- 
struction in our dwarf galaxy sample. It is also interest- 
ing to note that the increasing upper envelope in the age- 
mass distribution of Figure 2 (e.g., older clusters tend to 
be more massive) is a "size-of-sample" effect. This effect 
is evidence that significant m ass- independent destruc tion 
is not occurring in our data (Gieles & Bastian 2008). 



5.4. M, 



brightest 
V 



SFR 



We again use SLUG to test the effects of stochasticity 
on the results of the My" 6 cs — SFR diagram (Figure 
7). The same simulations as those in Section 5.1 are used, 
but this time we extract the brightest cluster formed in 
each iteration. The median brightest cluster is set as 
the simulation data point and the middle 68.3% of the 
simulations is set as the ler confidence interval. 

Figure 9 shows the <10 Myr and <100 Myr age ranges 
in the upper and lower panels, respectively. The open 
triangles are the SLUG simulated galaxies, and the red 
asterisks and blue x's are the <10 Myr and <100 Myr age 
range data points, respectively, with detected clusters 
above our completeness limit. The simulated M v " s 
values are set to a value of when the median of the 
simulations produced no clusters. The x-axis has been 
sorted by SFR and is increasing to the right. Simula- 
tions for the majority of the lowest SFR galaxies (left of 
UGC8833 and KKR3 for the <10 Myr and <100 Myr 
age ranges, respectively) produced very few detectable 
clusters in agreement with the lack of bright clusters ob- 
served in those galaxies. Data to higher SFRs in both age 
ranges show good agreement with the SLUG simulations. 

Due to the short time period of the <10 Myr age 
range there should be little effect from cluster destruc- 
tion, but there are five upper limit outliers that fall sig- 
nificantly below the expected relation in Figure 7. How- 
ever, these galaxies (UGC8508, NGC4163, NGC3741, 
UGC5139, and UGC8201) all fall within the expected 
scatter from the stochastic simulations. The agreement 
of these outliers with stochastic si mulatio n s sup ports 
the "size-of-sa mple" explanation of iLarsenl (J2002I ) and 
iBastianl J2008)'s conclusion that no outlier should signif- 
icantly fall below the trend. 

There are two galaxies that fall significantly above the 



M 



brightest 



SFR SLUG simulations: UGC5336 in the 



<10 Myr age range and UGC5692 in the <100 Myr age 
range. However, UGC5692 is shown not to be an out- 
lier if the <10 Myr age range data is used as discussed 
in Section 4.3. The <10 Myr outlier, UGC5336, shows 



moderate agreement with the M, 



brightest 



SFR trend, 



but shows poor agreement with the simulations. This 
cluster is a 98.2% outlier, i.e. 982 of 1000 SLUG simula- 
tions produced clusters fainter than the observed cluster. 



13 



4-4-4-4- 



<10 Myr 

A Simulated Data 
* Data With Clusters 
. Data W/0 clusters 



I I I I 



--cr)toco_K)cO"d-ir) 



- 1 A ± ± 



CO CM Kj «- 

O Cf. w ^ 

in co *- r- 

03 lO V "I 



i i i i 



<100 Myr 

A Simulated Data 

X Data With Clusters 

I Data W/0 Clusters 



jUj,4,UW T i 



Ss >. y 



-AAAAAAAAAAAAA 



to to 

to r-- 

- - CM (D ^- 



H 



x I 



i co c ,_ 7S ■ 

' " <3 CO D_ 
? 3 .O E O 



>CO^CO'tO>^lflKt 



) ^- cd lo 
) o co o . 

i n -1 r 



Increasing log SFR (M yr ') — > 

Fig. 9. — The simulated brightest cluster (My ls cs ) results 
from SLUG, where the top and bottom panels are the <10 Myr 
and <100 Myr age ranges, respectively. Each galaxy is presented 
on the ai-axis evenly spaced and sorted by SFR, with increasing 
values to the right. The open triangles are the SLUG simulated 
galaxies and the error bars are the 68.3% confidence inertvals of the 
1000 realizations of the simulations. The blue X's are the observed 
data and the downward arrows are the absolute V magnitude of 
the completeness limit for galaxies with no detectable clusters. The 
majority of the observed data from both age ranges fall within the 
expected scatter from the SLUG simulations, and therefore show 
good agreement. However, there are a few significant outliers above 
and below the simulations which might indicate a true variation in 
the cluster formation of these galalxies. 



5.5. Effects of Incompleteness 

Our completeness limit is only an approximate one and 
therefore it is possible that this limit can have an effect 
on the conclusions supported by our analysis. However, 
we will show in this section that the basic results are 
independent of the chosen completeness limit. 

The aspect of our analysis most affected by the com- 
pleteness limit is the upper limit data points. The upper 
limits represent galaxies without cluster detections and 
are computed directly from the completeness limit. The 
M v llg cs upper limits are computed as the complete- 
ness limit at the distance of the galaxies without cluster 
detections and the T upper limits are computed as the 
total mass of clusters below the completeness limit. 

The two main conclusions of our analysis are 1) true 
variations in the cluster formation of dwarf galaxies and 
2) that a mass-dependent cluster destruction model is 
dominant in our sample. The variation of cluster forma- 
tion in dwarf galaxies is supported by the significant out- 
liers in the F-TIsfr relation not accounted for by stochas- 
ticity. The supporting data are comprised of galaxies 
with cluster detections (e.g., UGC5336, UGC4459 and 



NGC2366) and no upper limits. The cluster destruction 
conclusions, on the other hand, are supported by some 
upper limits. However, only four of the eleven galax- 
ies used to draw these conclusions are upper limit data 
points. If the upper limits are ignored there would still 
be seven galaxies with cluster detections showing signif- 
icantly low total cluster mass all of which show evidence 
for the most massive clusters surviving. 

6. SUMMARY 

We have quantified the star formation and cluster for- 
mation properties of a large sample of dwarf galaxies. 
We use UBVRI data to derive the ages and masses of 
the cluster sample. Due to the low mass nature of the 
clusters in our sample we took extra care in the clus- 
ter parameter fitting using a Bayesian technique which 
specifically accounts for the stochastic sampling of stellar 
masses in low-mass clusters. 

The cluster sample was compared to previous studies 
by examining the T-Hsfr and M v " s os -SFR relation- 
ships in the age ranges <10 Myr and <100 Myr. This 
work provides the first examination of the these relation- 
ships in the low SFR regime. We found considerable 
scatter in our da ta compared to the T -Y^sfr. relation- 
ship proposed by iGoddard et al.1 (|2010f ). We also found 
moderate scatter in the M v " s es -SFR relation for our 

dwarf galaxies. In the My" s cs -SFR relationship some 
of our galaxies that lack clusters (upper limits) have pre- 
dicted cluster magnitudes 3-4 magnitudes brighter than 
our upper limits. In the T-T,sfr relationship NGC2366 
shows a significant deviation to lower T values and has 
als o been found to hav e unusually low cluster formation 
bv lBillett et al.l (|2002f ). Also, high T deviations are seen 
in UGC5336, UGC4459, UGC8760, and UGC5692. 

We model the expected scatter resulting from stochas- 
tically populating the CMF of our cluster sample using 
the SLUG software. We find that some of the outliers 
from the T-T,sfr and My 16 es -SFR relations are ac- 
counted for by the stochastic scatter, but there are galax- 
ies which show significant deviations beyond the stochas- 
tic scatter. The large outliers in the M v n& cs -SFR rela- 
tionship that fall 3-4 magnitudes below their predicted 
brightest clusters are within the stochastic scatter. Also, 
both UGC5692 and UGC8760 show significant scatter 
from the T-T,sfr relationship, but fall within the scatter 
seen in the SLUG simulations. NGC2366, on the other 
hand, shows significant deviations from the T-Y>sfr. re- 
lationship and the SLUG simulations in the <100 Myr 
age range. Furthermore, SLUG produced zero simulated 
cluster populations with T values as high as the observed 
data points for UGC5336 and UGC4459 in the <10 Myr 
age range. Inspection of the resolved star CMD of each 
galaxy's clusters and Ha images show correct age esti- 
mates. The cluster properties of these galaxies suggest 
that true variations in the clustered star formation are 
taking place in these low star formation rate galaxies. 

Although we find that the majority of individual galax- 
ies in the T-T,sfr relationship show moderate agree- 
ment (~ 2<r)with the stochastic modeling of SLUG in the 
<100 Myr age range, the galaxy sample as a whole shows 
systematic T values lower than the SLUG simulations. 
These lower values can be attributed to cluster destruc- 
tion effects, since this age range will be affected more 
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by cluster destruction. We constrain cluster destruction 
models by testing the mass-dependent model prediction 
that massive clusters preferentially survive cluster de- 
struction. If a mass-dependent model is dominant in our 
sample, the galaxies with low T values (i.e. low total clus- 
ter mass) should contain the brightest (i.e. most massive) 
cluster predicted by the M V " B es -SFR relationship. 

In our dwarf galaxy sample, 11 galaxies with signifi- 
cant scatter to lower V values in the <100 Myr age range 
all were found to have good to moderate agreement with 



the brightest cluster predicted by the M- 



brightcst 
V 



SFR re- 



lationship. In other words, galaxies exhibiting a dearth 
of total cluster mass in an age range where cluster de- 
struction becomes important all show evidence that the 
most massive clusters preferentially survive. This sug- 



gests that a mass-dependent cluster destruction model is 
dominant in our sample of dwarf galaxies. 
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TABLE 4 
Cluster Parameters 



Galaxy 


Cluster 


RA 


Dec 


U 


B 


V 


R 


I 


log(Age) 


log(Mass) 




(#) 


(J2000.0) 


(J2000.0) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(yr) 


(M Q ) 


ngc2366 


03 


07 28 40.4 


69 10 09.2 


21.01±0.04 


21.45±0.06 


21.01±0.02 


20.82±0.06 


20.28±0.03 


7.6 


3.4 


ngc2366 


05 


07 29 00.8 


69 14 08.2 


21.15±0.09 


21.78±0.10 


21.85±0.13 


21.82±0.26 


21.17±0.14 


7.6 


3.4 


ngc2366 


06 


07 29 15.3 


69 14 36.2 


20.44±0.03 


20.40±0.05 


19.97±0.01 


19.37±0.12 


18.71±0.01 


7.1 


3.7 


ngc2366 


18 


07 29 04.0 


69 13 01.3 


17.92±0.03 


19.13±0.05 


19.01±0.02 


19.11±0.12 


19.66±0.03 


6.5 


3.4 


ugc4305 


03 


08 18 54.5 


70 42 38.3 


20.17±0.11 


20.71±0.08 


20.63±0.06 


20.49±0.09 


20.18±0.06 


7.6 


3.8 


ugc4305 


23 


08 18 53.1 
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20.63±0.14 
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1.4 
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3.5 
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ic2574 
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10 28 44.6 


68 28 11.4 


17.87±0.05 
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18.41±0.06 
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ic2574 
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10 28 44.6 
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6.8 
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ic2574 
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10 28 37.3 


68 27 57.2 


18.69±0.04 


19.55±0.06 
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19.83±0.04 


7.2 
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ugc5692 


13 


10 30 35.0 


70 37 08.3 


19.93±0.11 


18.89±0.03 


18.46±0.03 


18.93±0.05 


18.64±0.04 


6.8 


3.1 


ugca292 


10 


12 38 40.0 


32 46 02.4 


19.65±0.04 


20.49±0.04 


20.34±0.04 


20.19±0.05 


20.39±0.03 


7.6 


3.9 


ugc8201 


08 


13 6 30.3 


67 41 51.1 


20.26±0.04 


20.60±0.04 


20.51±0.03 


20.30±0.05 


20.03±0.03 


7.8 


1.2 


ugc8201 


09 


13 6 22.6 


67 42 01.2 


19.49±0.03 


19.92±0.03 


19.90±0.02 


19.81±0.04 


19.52±0.03 


7.6 


1.1 


ugc8201 


10 


13 6 18.0 


67 42 12.8 


18.62±0.02 


19.08±0.02 


19.14±0.02 


18.79±0.03 


18.76±0.01 


7.4 


1.5 


ugc8201 


11 


13 6 17.6 


67 41 55.1 


21.94±0.12 


22.85±0.18 


22.50±0.21 


22.81±0.79 


22.09±0.30 


7.6 


3.4 


ugc8201 


12 


13 6 16.5 


67 42 04.8 


20.69±0.09 


21.04±0.08 
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7.2 


3.7 


ugc8201 


34 


13 6 29.9 


67 42 13.9 


20.73±0.05 


21.23±0.05 


21.22±0.05 


20.91±0.07 


20.81±0.06 


7.6 


3.9 


ugc8201 


38 


13 6 25.9 


67 42 03.7 


21.30±0.09 


21.09±0.06 


21.77±0.10 


21.07±0.10 


21.21±0.10 


7.1 


3.4 


ugc8760 


16 


13 50 51.7 


38 01 18.6 


19.15±0.02 


20.59±0.07 
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20.66±0.16 


6.5 


3.1 


ugc9128 


09 


14 15 55.0 


23 03 02.7 


20.44±0.06 


21.27±0.08 


21.27±0.12 


21.71±0.19 


21.58±0.24 


7.2 


3.1 


ugc9240 


12 


14 24 42.9 


44 31 42.4 


20.22±0.09 


21.26±0.19 


21.32±0.19 


20.65±0.24 


20.93±0.22 


7.6 


3.7 



Note. — The photometric properties of the cluster sample with an age less than 100 Myr.Wc report each cluster's identification number, 
right ascension, declination, UBVRI magnitude with error, the fitted age, and the fitted mass. 
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